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The behavior of hadronic matter at high baryon densities is studied within Ultra- 
J> 1 relativistic Quantum Molecular Dynamics (URQMD). Baryonic stopping is observed for 
Au+Au collisions from SIS up to SPS energies. The excitation function of flow shows 
strong sensitivities to the underlying equation of state (EOS), allowing for systematic 
studies of the EOS. Effects of a density dependent pole of the p-meson propagator on 
dilepton spectra are studied for different systems and centralities at CERN energies. 

OV 

*rH \ 1 Introduction 

> ■ 

The only possibility to probe excited nuclear matter in the laboratory are nucleus-nucleus 
reactions [Q. In particular when two heavy ions like Au or Pb collide most centrally, the 
combined system forms a zone of high (energy) density and high agitation of the involved 
constituents. The transient pressure at high density has specific dynamic implications, 
such as collective sideward flow. Hence, fundamental properties like the repulsion of 
the nuclear equation of state (EOS) are studied via event shape analysis of nucleons 
and clusters [QUIJ. The EOS at fixed temperature interpreted microscopically yields a 
density dependent potential modifiying the nucleon mass. At low densities this effect is 
similar as has been proposed by calculations on chiral limits of the Skyrme lagrangian||, 
the constituent quark mass|| and the chiral condensate (gg)0. Since (qq) should relate 
closely to hadron masses, the decay of short lived vector mesons, observed through the 
dilepton channel, is suggested as a promising experimental signal to investigate the gradual 
restoration of chiral symmetry. 

2 Ultrarelativistic Quantum Molecular Dynamics 

Since many important aspects of nuclear matter are not observable, numerical transport 
models are suited to test which assumptions are compatible to nature. The present model 
(URQMD) [|,[| includes explicitely 50 different baryon species (nucleon, delta, hyperon 
and their resonances up to masses of 2.11GeV) and 25 different meson species (including 
strange meson resonances), which are supplemented by all isosp in-projected states (see 
Table 1). Symmetries regarding time inversion, iso-spin, charge conjugation, etc. are 
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Table 1: List of implemented baryons, mesons and 
their resonances. In addition all charge conjugate 
and iso-spin projected states (and photons) are 
taken in and treated on the same footing. 
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Figure 1: a-i partial decay rates into 
specific channels. The average life- 
time is given by the inverse of the 
sum. Hence, in URQMD particles 
below resonance mass live longer, 
due to shrinking phase space. 



implemented in a general manner, e.g. all corresponding antiparticles are included and 
treated on the very same (charge-conjugate) footing. For excitations of higher masses 
a newly developed string model is invoked. It consistently allows for the population 
of all included hadrons from a decaying string. At low energies the dominant part of 
MM and MB interactions are modelled via s-channel reactions (formation and decays of 
resonances), whereas BB interactions are designed as exchange of charge, strangeness and 
four momentum in the t-channel. The real part of the baryon optical potential is modeled 
according to the Skyrme ansatz, including Yukawa and Coulomb forces. 

3 Creation of dense nuclear matter: stopping 

Baryonic stopping is a necessary condition for the creation of hot and dense nuclear 
matter. The key observable is the rapidity distribution of baryons. It is displayed in 
Figure |2] and |3| for heavy systems such as Au+Au and Pb+Pb at energies referring to 
three presently used heavy ion accelerators. In all cases gaussian rapidity distributions 
with peak around midrapidity y s ~ ±0.2 are found. However, the physical processes 
associated show characteristic differences: The average longitudinal momentum loss in 
the SIS energy regime is mainly due to the creation of transverse momentum whereas at 
AGS/SPS energies abundant particle production consumes a considerable amount of the 
incident beam energy. 

At CERN/SPS energies baryon stopping is influenced also by the formation time of 
strings which are excited in hard collisions. In URQMD baryons originating from a leading 
constituent (di-) quark at the string edges interact with (2/3)1/3 and mesons with 1/2 of 
their full cross sections during their formation time r. The sensitivity on this reduction 
is shown in Fig.|4] for the system S+S at 200v4GeV. The default calculation (including 
formation time) reproduces the data |TIJ fairly well whereas the calculation with zero 
formation time (dotted line) exhibits strongest stopping. A calculation with zero cross 
section within the formation time gives transparency. 
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Figure 2: Rapidity distributions for 
Au+Au collisions at SIS (lAGeV), AGS 
(10.6AGeV) and Pb+Pb at CERN/SPS 
energies (160v4GeV). All distributions have 
been scaled to the projectile rapidity in the 
center of mass frame. 
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Figure 3: Rapidity distributions for 
Pb+Pb at 160AGeV. The histograms la- 
bel from top to bottom: negative hadrons 
(h~), kaons (K~ , K + ), protons (b + — b~) 
and lambdas (A). The kaons are multiplied 
by five. 

In order to study this effect more closely the -y/s distribution for S+S reactions at 
SPS energies is shown in Fig.|5]. The collision spectrum exhibits two pronounced peaks 
dominated by BB collisions, one in the beam energy range and one in the low (thermal) 
energy range. Approximately 50% of the collisions, most of them around ~ 10 ± 5GeV, 
involve baryons during their formation time, whose cross sections are reduced by factors 
of 2/3 (referred to as di-quarks ) or 1/3 (referred to as quarks). This reduction leads to 
more transparency due to less collisions as compared to a calculation without reduction. 
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Figure 4: Sensitivity of the rapidity distri- 
bution for S(200v4GeV)S on the constituent 
(di-) quark cross section. 
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Figure 5: ^ distributions for baryon 
baryon collisions in central reactions of 
S+S (right) at SPS energies. 
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Figure 6: Excitation function of the total 
directed transverse momentum transfer px- 
dir for Au+Au. URQMD calculations in- 
cluding a hard EOS (full squares) are com- 
pared to the predictions of cascade calcu- 
lations (full circles). 
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Figure 7: Mean directed transverse mo- 
mentum as a function of the scaled rapid- 
ity. The transverse flow < p x > (y/y pro j.) 
scales with the mean transverse momentum 
< p t >, i.e. directivity does not depend on 
the bombarding energy. 



4 Probing the repulsion of the EOS: flow 

The creation of transverse flow is strongly correlated to the underlying EOS [|TJ. In par- 
ticular it is believed that secondary minima as well as the quark-hadron phase transition 
lead to a weakening of the collective sideward flow. The occurrence of a phase transition 
should therefore be observable through abnormal behavior (e.g. jumps) of the strength 
of collective motion of the matter ]T2[] . Note that URQMD in its present form does not 
include any phase transition explicitly. In Fig.^| the averaged in plane transverse momen- 
tum is displayed for Au+Au from 0.1 to 4 AGeV incident kinetic energy. Calculations 
employing a hard EOS (full squares) are compared to cascade simulations (full circles). 
In the latter case only a slight energy dependence is observed. In contrast, the calcula- 
tion with a hard EOS shows strong sensitivity. Here, the integrated directed transverse 
momentum per nucleon is more than twice as high as for the cascade calculation. This 
indicates the importance of a non-trivial equation of state of hadronic matter. 

The amount of directed transverse momentum scales in the very same way as the total 
transverse momentum produced in the course of the reaction. Hence, the directivity de- 
pends only on the reaction geometry but not on the incident energy. This is demonstrated 
through Fig. [7|, where the mean p x as a function of the rapidity divided by the average 
transverse momentum of all particles is plotted. 

5 Temperature dependence of the EOS: photons 

Semiclassical cascade models in terms of scattering hadrons have proven to be rather 
accurate in explaining experimental data. Therefore it is of fundamental interest to extract 
the equation of state from such a microscopic model, i.e. to investigate the equilibrium 
limits and bulk properties, which are not an explicit input to the non-equilibrium transport 
approach with its complicated collision term (unlike e.g. in hydrodynamics [|T2| , p~3|] ) . In 
Fig.|8]the thermodynamic properties of infinite nuclear matter are studied within URQMD. 
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Figure 8: 'EOS' of infinite nuclear matter 
function of the energy density versus 
temperature at fixed net-baryon density of 
p B = 0.16/fm 3 in URQMD (symbols). The 
curves refer to analytical forms of the EOS, 
i.e. a Hagedorn-gas (top), a quark-gluon 
plasma (middle), and an ideal gas of nu- 
cleons and ultrarelativistic pions (bottom). 
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Figure 9: Transverse momentum spectrum 
of directly produced photons in Pb+Pb 
collisions at 160v4GeV calculated with 
URQMD. The resulting spectrum is com- 
pared with hydrodynamical calculations. 
In all models the processes ni] \— > rcy, 
up i— > 7T7 and tttt i— > p7 are considered as 
photon sources. 



Infinite hadronic matter is simulated in URQMD by constructing a box of 250fm 3 
volume with periodic boundary conditions. According to saturation density, nucleons are 
initialized randomly in phase space, such that a given energy density is reproduced. After 
the system has equilibrated according to the simulation with URQMD the temperature 
is extracted by fitting the particles' momentum spectra. Alternatively, the temperature 
can be extracted from the relative abundances of different hadrons, e.g. the A/N ratio. 

In Fig.p the result of this procedure is compared to various analytic forms of the EOS. 
While the EOS of a Hagedorn gas and a QGP yields energy densities e ~ lGeV/fm 3 at 
T = 150MeV the temperature dependence is much smaller in URQMD. It yields about 4-5 
times less energy density, being in fair agreement with a gas composed of nonrelativistic 
nucleons and ultrarelativistic pions. It remains to be seem whether a reparametrization of 
the resonance continuum in the Hagedorn model as suggested in Ref. would resolve the 
deviation as compared to URQMD. On the other hand, beyond T ~ 200MeV the energy 
density rises much faster than T 4 approaching even the QGP value of e ~ 10GeV/fm 3 
around T = 300MeV. This indicates an increase in the number of degrees of freedom. It 
may be interpreted as a consequence of the numerous high mass resonances and string 
excitations, which seem to release constituent quark degrees of freedom (but, of course, no 
free current quarks as in an ideal QGP). Investigations of equilibration times and relative 
particle and cluster abundances are in progress. Moreover, the admittedly poor statistics 
have to be improved, in order to study the high temperature behavior. 



Experimentally, the EOS can be accessed by measuring electromagnetic radiation [14]. 
In Fig.RI the direct photon production from meson+meson collisions in Pb+Pb collisions 
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Figure 10: Dilepton mass spectrum for 
p+Be at 450GeV. The calculation includes 
Dalitz decays and conversion of vector 
mesons (see also legend for S+Au). The 
sum of all contributions (solid curve) is 
folded with the CERES mass resolution. 
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Figure 11: Dilepton mass spectrum for 
S+Au at 200AGeV (see also legend for 
p+Be). Here no in- medium modifications 
of the p propagator is considered. Around 
M ~ 400MeV two points are missed by 
about two standard deviations. 



at 160AGeV is shown. Here, only mesons stemming from string decays are included. 
Elastic meson+meson scattering with a e \ = 15mb (independent of a/s) was considered. 
The result is compared to calculations within the 3-fluid model scaling and Landau 
expansion with Tj = 300MeV. 

6 In medium masses: dileptons 

In Fig.|10| and [Ll] calculations of dilepton spectra with URQMD are show for p+Be and 
S+Au. Dilepton sources considered here are Dalitz decays (tt°, rj and u) and vector 
meson decays (p, uj and 0). Dalitz decays of heavier meson and baryon resonances are 
included explicitely via their emission of p mesons (assuming vector meson dominance). 
In order to avoid double counting, the p mesons from 77's, and a/s are excluded from the p 
contribution. Pion annihilation is included dynamically into the contribution of decaying 
p mesons (7r + 7r~ \— > p \— > e + e~). 

While the result for p+Be agrees well with the published data from CERES/SPS 
|T7|1 , two points around M ~ 400MeV are missed by about 2a for S+Au. Speculations 
about the origin of this deviation propose electromagnetic bremsstrahlung, annihilations 
of pions and a modification of the p meson propagator due to a gradual restoration of the 
chiral symmetry. The contribution of pion annihilation to the p-peak (7r + 7r~ h- > p) is only 
40% for S+Au. Major additional sources are decays of heavy baryons A*/A r * 1— > Np as 
proposed in Ref.JT^] and meson resonances (see also Fig.|TJ) : 
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In Ref. [^,[7| a linear dependence of the p°/u; pole position as a function of the nuclear 
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Figure 12: Dilepton mass spectrum for 
S+Au at 200v4GeV. The curves label sim- 
ulations with pole shift according to the 
creation density (top), the decay density 
(middle) and without pole shift (bottom). 
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Figure 13: Dilepton mass spectra for 
Pb+Au at 200AGeV. The curves label cal- 
culations for central (cent.) and inclusive 
reactions (incl.) with (mod.) and without 
a pole shift (free). 



density p has been suggested: m p o(g/g ) = m p o(0) (1 — Xg/go)- Here p denotes the 
ground state density of nuclear matter, and A = 0.18, in agreement with various other 
calculations. Since the restriction to low densities may not be suitable for heavy ion 
collisions, the following extrapolation towards higher densities is taken: 



fn p o(g/g ) 



m p o (0) 
1 + A g/g 



(2) 



In Fig.|12] an application of Eq.(|2|) is made to calculate a dielectron mass spectrum for a 
density dependent vector meson pole. In URQMD the p meson pole position is shifted 
according to the density at which the p meson decays, i.e. eventually converts into e + e~ 
(middle curve). Note that this procedure is equivalent to a "shining" description where the 
p constantly emits e + e~ pairs according to the rate dN ee /dt = T(p i— ► e + e~). This result 
yields only a small enhancement around M ~ 500MeV as compared to the calculation 
without pole shift (bottom curve). On the other hand, the data can nicely be reproduced, 
if the strong (unphysical) assumption is made, that the pole at the decay point (p i— > e + e~) 
is shifted according to the creation density (upper curve). This would be a neglection of the 
finite decay length. The discrepancy of a calculation without decay length (Ad cc = Ofm) as 
compared to the result including the decay length is driven by two reasons: i) The increase 
of the p lifetime (~ 7fm/c) below its resonance mass in the region M ~ 0.3 — 0.5GeV 
(where a dilepton excess in S+Au is reported |17|]) lowers the decay density down to 
(g) ~ 0.2^ o for S+Au (or 0.3f? o for Pb+Au). ii) An enhancement of the decay length leads 
to an increase of reabsorption. Hence, the radiation path for p i— > e + e~ is substantially 
truncated. 

The result for Pb+Au is shown in Fig.|T3" 



are in fair agreement to the preliminary observation from CERES [19 



Both calculations for inclusive reactions 

For central events 



an enhancement due to a nonlinear in-medium p-decay contribution is predicted. The p 



bump is located either at ~ 500MeV with pole shift or at ~ 700MeV without modification. 
Note that the deviations induced by the pole shifts are for both systems and centralities 
lower than the statistical errors. Hence, the interpretation of the experimental data gives 
the following impression: The data for light systems such as p+Be and p+Au (see also 
Ref . pUf ) as well the data for heavy systems both for inclusive and central reactions are 
reproducible whithout mass shifts. In contrast the central data for S+Au exceed the 
URQMD calculation around M ~ 0.4GeV by about two standard deviations 

7 Summary 

Studies of the equation of state and consequences of gradual restoration of the chiral sym- 
metry are presented using a microscopic phase space model including 75 hadron species 
and strings. The directed transverse momentum shows strong sensitivities on the under- 
lying EOS: There is only a small increase in the cascade calculation, whereas it scales 
linearly with the average transverse momentum for a hard equation of state. Hence, mea- 
surement of the excitation function of the transverse directed flow allows for a systematic 
study of the EOS. 

The increase of the p lifetime at masses M ~ 0.4GeV where a dilepton excess in S+Au 
is reported |T7| lowers its average e + e~ emission density substantially. Therefore it seems 



that the density dependence of the p pole alone does not suffice to explain the dilepton 
excess in S+Au. Hence, other effects - such as the temperature dependence of (qq) or 
additional sources - might be required. 
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